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ABSTRACT 


An analysis is presented of the thrust requirements to compensate for 
gravity gradient torques and aerodynamic drag and to accomplish annual 
precession of rotating manned space stations. Analysis and design studies 
of electrically heated and solar -heated reaction control systems which use 
carbon dioxide as the working fluid are reported. Experimental data obtained 
by transient and steady state performance testing of an Ohmjet electrically- 
heated pulse rocket are discussed. Conclusions and recommendations are 
presented regarding the applicability of Ohmjet and Heliojet thrust systems 
for reaction control of manned space stations. 



I INTRODUCTION 


The objective of this program has been to conduct analytical and design studies 
to aid in the selection and design of precession control and orbit maintenance systems 
for a rotating space station, utilizing electrically-heated and solar -heated pulse rockets. 
Also experiments were conducted to determine the performance parameters of an Ohmjet 
electrical resistance -heated rocket of a type applicable for precession control. 

The effort was divided into the following five tasks: 

(a) Analysis and design studies to determine the required performance parameters 
for thrust producing systems. 

(b) Analysis of the design parameters for an Ohmjet thrust system. 

(c) Analysis of the design parameters for a solar -heated (Heliojet) thrust system. 

(d) Experimental performance evaluation tests of a prototype Ohmjet thrust unit 
which had been previously designed and fabricated by Exotech Incorporated. 

(e) Development and improvement of the Ohmjet thrust unit. 

The general reasons for interest in the application of electrically-heated and solar - 
heated pulse rockets for attitude control and orbit maintenance of manned space stations 
are: 

1 . The use of waste carbon dioxide from a manned spacecraft for reaction control 
would reduce or eliminate the requirement for additional chemical propellants. 
However, the availability of CC^ for thrust production would be contingent on a 
decision not to regenerate CC^ into oxygen because of system considerations 
such as short mission duration or lack of reliable oxygen regeneration equip- 
ment. 

2. The disturbing torques and forces on a manned space station include a fairly 
steady "base load" requirement to compensate for gravity gradient torques 
and aerodynamic drag and to provide precession at a steady rate of one revolu- 
tion per year for sun -orientation. These base load requirements can poten- 
tially be met by using waste carbon dioxide which would be generated by the 
crew at a fairly constant rate. Peak loads on the attitude control system which 
would be imposed by docking forces, motion of the crew, etc. could be met by 
a separate system incorporating chemical rockets and flywheels or gyros. 
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3 . The concept of using electrical or solar energy to heat a working fluid for re- 
action propulsion has broad potential for application to spacecraft if waste 
fluids can be used effectively or if a low molecular weight propellant such as 
hydrogen or ammonia can be used, so that the specific impulse may exceed 
that attainable by chemical rockets. The technology of such devices needs 
further development. 

It is hoped that the analysis, design and experimental efforts reported herein will 
aid in defining the directions for future development of the technology of reaction control 
systems using pulsed operation of electrically-heated and solar -heated thrust devices. 


2 



II ANALYSIS OF DISTURBING TORQUES AND THRUST REQUIREMENTS 


A . Objectives and Basic Assumptions 

The purpose of this analysis is to provide a basis for the preliminary design of 
the reaction control system for a rotating sun-oriented space station . The function of 
the reaction control system is to point the spin-axis toward the sun and to compensate 


for aerodynamic drag in order to maintain orbit altitude . 

Wherever possible the analysis is made general enough to render it applicable 
over a range of orbit and station parameters . For the purpose of making numerical 
estimates, however, the following specific values are used: 

Circular Orbit Altitude 300 nautical miles 

Orbit Inclination 30 degrees 


Spin Axis Moment of Inertia 
Transverse Moments of Inertia 


7 2 
I z = 1-5x10 slug-ft 

I x = l y = 1 -05 x 10 7 slug-ft 2 


Station radius (maximum moment arm 
for thrust application) 


75 feet 


Spin rate ° lvrM 

The reaction control system must have sufficient capacity to counter all external 
disturbances which alter station momentum (both linear and angular). These disturbances 
can be categorized as follows: 

I . Orbit Maintenance 

(a) Aerodynamic drag 

(b) Solar radiation-pressure (negligible for space stations) 

II . Attitude Control (Sun-Orientation) 

(a) Yearly precession around the sun 

(b) Gravity- Gradient 

(c) Aerodynamic torques 

(d) Solar-radiation pressure torques 

(e) Interaction with earth’s magnetic field 

(f) Meteoroid impacts 


The relatively low altitude and the fact that this is a spinning station makes 
some of the above disturbances dominant in comparison to the other . It was thus found 
that the largest single source of disturbance is that due to gravity gradient torques; the 
second largest disturbance source is that due to the yearly motion Of the earth . 
around the sun. The amount of average thrust required for orbit control purposes 
was found to be comparable to the average thrust needed to accomplish annual 
precession . 

The major disturbing torque produced by gravity gradient and the torque required 
to precess the station spin axis by one revolution per year (the annual precession 
torque) are both amenable to detailed analysis . In both cases , orbit and station 
design parameters entering into the analysis are known and remain relatively constant 
and the analysis can be based on well established laws of mechanics . Therefore the 
analysis can be rigorous and the results can have general utility in defining the pre- 
dominant long term (secular) torque and thrust requirements . The analysis was 
intentionally conducted in this general manner in order to estimate "base load" re- 
quirements and to avoid consideration of short time transient disturbing effects such 
as docking loads , motion of crew, etc . which can be separately counteracted and 
which involve detailed system and operational assumptions . 

B . Comparison of the Relative Magnitude of Disturbing Torques and Thrust Requirements 
It is convenient to use the annual precession torque as a reference and to compare 
all other torques with it . The torques are expressed in terms of the average thrust 
required at a moment arm of 75 feet from the center of mass . 

(a) Annual Precession Torque 

Motion of the earth and therefore the space station around the sun 
introduces the requirement that the station axis of symmetry, which is also 
the spin axis , must precess in inertial space at a rate of one revolution per 
year in order to maintain sun orientation which would permit use of solar 
power and simplify thermal control . The magnitude of the average torque 
needed to produce this annual precession is given by the gyroscope equation: 
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= spin rate about sun-pointing axis 
cO - 1 revolution per year 

For the numerical values previously defined 
T p = 0 .94 lb -ft 

and the average thrust at a moment arm of 75 ft is 
F p =0.0125 lbs 

The utility of the precession torque as a reference stems from the 
fact that it is independent of orbit parameters , i.e. altitude, or inclination, 
and that it is constant with time . 

It is noted that the annual precession torque vector must lie in the 
ecliptic plane and be perpendicular to the earth-sun line in order to produce 
precession around an axis perpendicular to the ecliptic plane . Therefore 
thrust must be applied toward or away from the sun at points above or below 
the ecliptic plane passing through the station center of mass . Since the 
station is rotating around the earth-sun line, it is impossible to produce a 
steady annual precession torque by continuous operation of a thrust device 
which is fixed to the space station . Therefore the thrust must actually be 
applied in pulses while the thrust device is above or below the ecliptic plane . 
(b) Gravity Gradient Torques 

We consider next the magnitude of gravity gradient torques . As derived 
in the course of this study, the average value over any one orbit of the secular 
component of the gravity gradient torque is given by the expression 

t g <W D s < 2 > 

where s orbital period 

D_ is a function of orbital parameters and its value depends upon initial 

u 

launch conditions . In general , Dg varies between zero and 1 with a period of 
about 54 days . However, for the assumed orbit inclination of 30° it is always 
larger than zero . For the present purpose of comparing magnitudes we will 
assume an average value of 0 .6 . 


As shown in reference 1 , for a spinning body in which , 

stability of rotation requires that 1^ be about 1 .5 to 2 times larger than L^. . 

Wc assume 

Z 

Hence 


>- =2, x 


- Sc 


( 3 ) 


l Z *X ' z 

To compare gravity-gradient torques with the precession torque we 
evaluate their ratio: 


t g 

"tT 'o 


D, 


( 4 ) 


For P in minutes and in RPM , equation (4) becomes 
0 z 


-( 


2 x io 5 ) D : 


(5) 


p ^Y> 

o AX z 


Evidently, the relative magnitude of these torques is influenced only by 
the spin rate and orbital period . For the particular case considered here , P Q 
varies only slightly with altitude, i .e . , it is 95 minutes at 300 n .miles and 
92 minutes at 200 n .miles . Using -T> z = 3 RPM we get a maximum ratio 

(for D s 1) of: 


(*) 


= 7.27 


max 


and an average ratio (for Dg - 0 .6) of: 


\ P / avc 


= 4.36 

, _ average 

Thus , the secular gravity gradient torque is substantially larger than the 

torQue required for suniuil precession • 

It is worth noting that the above will be true for all practical configurations 

f a rotating sun-oriented station where the function of rotation is to provide artificial 
cavity. For example, as a limiting case,assume that it is desired to provide a full 
g of centipetal acceleration at a radius of only 32 feet, fliis would require a spin rate 
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of 10 RPM and for die same range of orbital altitudes the average value of 
the secular gravity gradient torques would be 1 .32 times as large as the 
orbital precession torque. This ratio forD^ = 1 would be 2.2. Hence the 
validity of the basic conclusion is not significantly affected by changes in 
station design parameters . 


(c) Orbit Control 

We consider next the average thrust required to maintain orbit 
altitude . Although for completeness , radiation-pressure was listed as a source 
of disturbance, it is of no practical relevance because in the altitude range 
considered it is at least two orders of magnitude smaller than aerodynamic 


pressure . The principal cause of orbit decay is therefore air drag . 

It is assumed that the drag force is directed along the negative velocity 
vector and equals 

F D = M C DV >(h)V (6> 

A s projected area normal to velocity vector 


p(h) s air density as a function of altitude 
V - satellite velocity 


We will assume free molecular flow and completely diffuse reflection of 


molecules rebounding from the surface . Under these conditions , and neglecting 
the thermal velocity of impinging molecules by comparison with the orbital station 
velocity, the drag coefficient is: 

C D = 2 


To evaluate the projected area, we assume the station to be an annular 
disk of inner radius r and outer radius with its axis pointing toward the sun. 

Denoting by A the area of the ring: 

A = TT(t 2 2 - r t 2 ) (7) 

The projected area A p in the direction of the velocity vector varies from 
zero to A over each quarter orbit because of the variation of the angle ^ between 
the sun line and the velocity vector . For a circular orbit in the ecliptic plane. 

A = A cos <P 


7 



( 9 ) 


The average value of for a nearly circular orbit is: 
Ap A = f A /cos,* | A 


The average value of the drag force therefore becomes: 

F D = 2(r 2 2 - r l 2> / >(h > ^ 


( 10 ) 


The largest uncertainty in estimating the drag force is due to the lack 
of definitive data on the daily average of air density at the altitudes of interest . 
As shown in reference (2) there is a spread of one order of magnitude in the air 
density values reported by various investigators for any one altitude . We use 
data based on reference (3), viz . 

at 200 n .miles 
at 300 n .mil es 


-14 , 3 

— 3 x 10 gm/cm 

•15 3 

/o 3 x 10 gm/cm 


To estimate , we use r 2 =75 ft, ^=65 ft. The average drag force is 


then found to be: 

h = 200 n .miles 


D = 0.083 lbs 


h = 300 n .miles 


F D = 0 .008 lbs 


The average thrust to provide annual precession torques was previously 
found to be 0 .0125 lbs . At an altitude of 300 n .miles the average thrust required 
to compensate for drag is therefore approximately equal to that needed to provide 
annual precession torque . However at 200 n .miles the drag effect is about 6 .6 
times larger and is therefore comparable to the peak values of gravity gradient 
torque , with respect to thrust requirement . 

(d) Aerodynamic Torques 

The preceding calculations of air drag can be used to estimate aerodynamic 
torques on the station . These would be due to a displacement Ar of the center of 
pressure from the center of mass and could be calculated from 


T = F 
D D 


Ar (11) 

The value of Ar depends upon the size and geometry of the station and 
might be expected to be a few percent of the maximum radius , say: 

Ar = 0.02r 

max 
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Hence 


T D =0 - 02F D - r max 


(13) 


Using equation ( 13) the comparison between annual precession torque 
and aerodynamic torque is as follows: 


h = 300 n .miles 


1EL_ =0.013 


h = 200 n .miles 


D =0.133 


In either case , therefore , aerodynamic torques are substantially 
smaller than the required annual precession torque. 


(e) Solar Radiation -p res sure 

As previously noted, solar radiation pressure is a few orders of 
magnitude smaller than aerodynamic pressure at the altitudes of interest. 

Hence , the effect becomes negligible by comparison with the other disturbing 
torques . 

(f) Magnetic Interaction Torque 

Since the magnitude of torque due to interaction with the earth's 
magnetic field depends upon the amount of magnetic materials in the station 
and the size and location of current loops , it can be assumed that precaution 
in design will keep these torquesi down to a low value. In general, these torques 
may be expected to be less than the aerodynamic torques , unless the station were 
specifically designed to use magnetic torques for control . 

(g) Meteoroid Impacts 

Lack of definitive data on meteoroid fluxes makes a realistic estimate 
difficult. However, based on prior analyses, it is believed that the torque require- 
ments to compensate for meteoroid impacts will be negligible . 

(h) Summary 

The data developed above can be summarized in tabular form by showing 
the various average thrust levels as a fraction of the thrust needed to provide 
for annual precession . 
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Disturbance 

Fraction of Annual Precession Thrust 
h = 200 n .miles h = 300 n .miles 

Annual Precession 

1 

1 

Aerodynamic Drag Force 

6.6 

0.64 

Gravity Gradient Torque 

4.5 

4.36 

Aerodynamic Torque 

0.13 

0.01 

Magnetic Interaction Torque *1 — 

F - 

Meteoroid Impact I 

0.2 

0.2 

R aviation -Pressure Torque 3 — 

TOTALS 

12.43 

6.21 


The totals shown above indicate that propellant requirements would be 
twice as large at 200 n .miles as compared to an orbital altitude of 300 n .miles . 

It is to be noted, however, that since the difference is caused primarily by the 
increased drag on the station , and since this estimate is based upon assumed 
values of density which are not very reliable , the factor of two indicated above is 
rather uncertain . Even though there is no agreement on the absolute value of 
air density at any one altitude , it is probable that the average density changes by 
about one order of magnitude between 200 and 300 n. miles . Hence, from the 
point of view of minimizing propellant weight for reaction control purposes , 
the orbital altitude should be 300 n .miles . Assuming that this would be the 
design choice, the largest single source of torque would be gravity gradient. 
Furthermore , gravity gradient and annual precession torques would account for 
at least 87% of the total propellant weight required, even if no attempt were made 
to use a tt itude control impulses to accomplish drag compensation at the same time . 
C . Gravity Gradient Torques 

A preliminary analysis of gravity gradient torques on a rotating, sun-oriented 
station is reported in Topical Report TR-001, (Ref. 4). The following is intended to 
summarize Ref. 4 and the additional analyses performed during the program . 
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It is estimated that a saving of about 20% in propellant consumption can be 
achieved by applying the control impulses with a single nozzle at a time when station 
rotation has brought this nozzle to an angular position in the plane of the station which 
is 90 degrees away from the axis about which the torque is to act . From the point of 
view of evaluating disturbance torques , this implies that only the total torque magnitude 


is relevant, regardless of where this vector might lie in the plane of the station . The 
analysis of gravity gradient torques was therefore made in terms of the total torque 
magnitude , rather than separate components in the plane of , and perpendicular to , 
the ecliptic plane . 

The instantaneous value of the total gravity gradient torque is given by 

(14) 


J T = I (I z' I X )Sin2l % 


GM = gravitational constant for the earth 
Rq = orbital radius 

I^ = moment of inertia about spin axis 


^ 1 2 moment of inertia about axes in the plane of the station 

~ between the sun -pointing spin -axis and earth's radius vector . 


Since ^ varies with the instantaneous value of the various orbital parameters , 

the total gravity-gradient torque will be time varying . However , because the period 

of the satellite in its orbit around the earth is much shorter than the period of the 

orbit around the sun , the gravity gradient torque can be defined as a product of two 

terms , one of which is constant over any one orbit . Integrating the torque over a 

complete orbit, the following expression was derived for the angular momentum per 

orbit due to the total gravity gradient torque: 

12tt 

h t = P — <I Z -Ix> f (D) (per orbit) ( 15 ) 

where 

f(D)=D*$/2)(l-D 2 ) 

D =y - sin 2 i sin 2 ^ 

i = angle of inclination of the orbit plane to the ecliptic plane 
at launching and which changes slowly as orbit regression 
occurs . 

- angle between the spin axis and the intersection of the orbit 
plane with the ecliptic plane , which also varies with time . 
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The parameter f (D) is a function of the various orbital parameters and changes 
slowly with time . Its definition requires a knowledge of orbit inclination , altitude 
and time of launch . If these parameters are specified, the propellant consumption 
required for torque compensation , or the station angular misalignments can be 
evaluated as a function of time . 

For the assumed orbital parameters , the maximum value of f (D) is 1.2. Although 
theoretically f (D) can also reach zero (when i and^Fare 90°), this could not be the case 
for the assumed value of i = 30° . 

A similar approach was taken to evaluate only the secular component of the total 

gravity gradient torque . Angular momentum over one orbit is in this case given by 

~ 2 

H g — <I Z - l x ) Dg (per orbit) (16) 

where 

Dg = 2 sin i siny^||l - sin^ i sin^|P 

D c is also a function of the time varying orbital parameters , but differs from 
f(D) . The maximum value of Dg is 1 and , for the assumed value of i = 30 , will 
always be larger than zero . 

D . Station Angular Misalignment 

Assuming that only gravity-gradient torques are acting on the station , it is of 
interest to estimate the rate at which the sun-pointing axis would be misaligned if no 
control torques were applied . It is also relevant to evaluate the separate effects of 
periodic and secular gravity-gradient torques . 

Since the analytical relationships do not lend themselves to a separate evaluation 
of the periodic component of the total torque , an exaggerated case is used to obtain an 
upper limit. It was thus assumed that the magnitude of the periodic component is equal 
to the magnitude of the total torque with f(D) at the maximum value of 1.2. The periodic 
component acts at twice the orbital frequency and the peak value of angular misalignment 
was found- to be 0.115 degrees during one quarter of an orbit. Hence, for the assumed 
numerical values , the periodic torques will not greatly influence the attitude control 
problem . 
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Angular misalignment due to the secular torques , with D g at the maximum value 
of 1 , was found to be 


% 


max 


= 0.3 degrees per orbit 


Figure II- 1 shows the accumulation of angular misalignment for the above maximum 
rate conditions . It is therefore seen that the accumulated error over a day would be 
less than 5 degrees even at the worst conditions . The superimposed oscillation due 
to periodic torques is also seen to be negligible by comparison . 

Figure II-2 illustrates the effect of spin rate upon angular deviations . It is of 
interest to note the significant errors which would have been produced by the same 
gravity gradient torques if the vehicle were not spinning, name 13 .5 degrees per 
quarter orbit (22 .5 minutes) as compared to 0 .075 degrees in the same period for a 
spin rate of 3 RPM , In particular , if the station were not spinning , additional control 
means , e .g . , a momentum wheel would be required to counter the periodic torques 
since they would otherwise cause appreciable misalignments . 

The angular stiffness of the station , illustrated by the above numerical estimates , 
indicates considerable latitude in the programming of control impulses . To illustrate 
this further , it was assumed that control impulses are applied continuously during the 
time when secular gravity gradient torques are near their peak value (D g = 1) but that 
the level of the applied thrust is less than that required to completely counter the 
gravity gradient torque . Assuming the control thrust to be 90% of the peak gravity- 
gradient torque and assuming also that the latter varies sinusoidally around the peak 
value , it was found that the net misalignment during the time when the control torque 
is less than the applied torque would be 2 .6 degrees . This error would be accumulated 
over a period of 7 .8 days . This calculation is based ijpon a computed value for the 
long-term period of gravity-gradient torques of 54 days . 

E . Estimates of Propellant Weight 

This estimate assumes that a propellant having a specific impulse of 300 seconds 
is carried aboard the station to provide reaction-jet control for station orientation and 
orbit control . For the particular station parameters as used here , and assuming an 
orbital altitude of 300 n .miles , the required weight of propellant for a one year period 
is estimated as follows: 
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FIGURE II 
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(a) It is assumed that angular orientation requirements will permit a mis- 
alignment of 0 .2 degree or more . Hence , no control is provided for the periodic 
component of gravity-gradient torques . 

(b) Angular momentum accumulation due to the secular gravity gradient 
torques depends upon the time variation of Dg. Since Dgi 1 , an upper limit is 
obtained by letting D = 1 . Propellant weight due to secular gravity gradient 
torques is thus found to be 

W c < 6 ,480 lbs per year 

Since D cannot be zero , the average value of D over long periods of 

s ^ 

time might be expected to be 0 .6 or more . Hence 
3 ,880 < Wg < 6 ,480 lbs/year 

(c) The amount of propellant to produce only the annual precession rate is 
calculated to be 


Wp = 1 , 350 lbs per year 

The annual precession torque required is unidirectional and cannot , 
therefore , be added directly to the secular gravity-gradient torque . For , if 
during some portion of the year the annual precession torque adds to the gravity 
gradient torque , during other times it will tend to compensate for gravity gradient . 

(d) The combined effect of drag causing orbit decay , and the various other 
disturbance torques is to require additional propellant weight comparable to 
that required for annual precession torque , i .e . 

^W. , 9? 1,350 lbs per year 

Based on the above it is estimated that chemical propellant weight for the 
300 n .mile altitude would be about 6 , 000 lbs per year . If a chemical bipropellant 
rocket system were used, its total weight including propellant tankage and thrust 
devices would be approximately 7500 lbs for a one year life with no system 
redundancy . 

Assuming a space station weight of 120,000 lbs, a chemical reaction 
control system would require about 6 % of the station gross weight for a 300 nautical 
mile orbit . For lower altitudes , this fraction would be higher . 
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F . Methods of Thrust Application 

The manner in which control is to be applied to a large , rotating station should be 

governed by the following considerations: 

(a) For a 300 n .mile altitude , about 87% of the thrust demand is predictable 
as regards its variation with time . This includes the requirement to correct 
for secular gravity gradient torques and to apply annual precession torque . 

(b) The angular stiffness of a rotating station provides a great deal of 
latitude in the programming of applied thrust. Thus , if an analytical prediction 
of thrust demand is used , it need not be extremely accurate . More significantly , 
however , station attitude control would not be dependent upon precise control of 
thrust magnitude or duration . 

(c) Station rotation can be used to permit firing of a single nozzle at the 
desired angular position. As previously described, this results in a reduction 
of propellant consumption . 

(d) Thrust requirement for orbit maintenance purposes is about 15% of the 
total for a 300 n .mile altitude . Since a single nozzle does not produce a pure 
couple on the station, it is possible to program the attitude control impulses 
from a single thrust device so that when the station is moving away from the 

sun , the impulses occur with a different frequency from when it is moving towards 
the sun . This will produce a net force on the station, opposing the drag force 
which causes orbit decay. If two-thrust devices are used, one on the sunward 
side and one on the shadow side , the air drag compensation can be accomplished 
even more efficiently as a byproduct of attitude control impulses . 

Since the angular position (in the plane of the station) at which thrust is 
applied is one of the control parameters , thrust duration must be sufficiently short 
so as not to spread it over too large an angle by virtue of station rotation • Other- 
wise the effective moment arm of the impulse would be reduced. Figure II -3 shows 
the efficiency of thrust application as a function of the ratio of firing time to the 
period of station rotation . For the assumed rate of 3 RPM and for a firing time 
of 2 sec the efficiency is seen to be about 92%. 
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Effic ie ncy 



Ill DESIGN STUDIES OF AN ELECTRICALLY-HEATED 


PULSE ROCKET (O HMTET) THRUST SYSTEM 


A. General Discussion - 

The inherent advantage of using waste carbon dioxide from a manned space station 
as the propellant for attitude control can only be realized through an efficient means of re- 
covery and utilization of this waste product. It is assumed that a molecular sieve bed is 
used to absorb C0 2 - Two basic methods for recovery of C0 2 are described and analyzed 

in this section; they are: 

Method I: Removing the carbon dioxide from the molecular sieve bed by heating 
the bed until the CC> 2 partial pressure is raised to 14.7 psia; compressing this gas to ten 
atmospheres and storing it in a storage bottle, then exhausting it periodically, as required, 

through a small Ohmjet. 

Method II: Removing the carbon dioxide from the molecular sieve bed as above, 
but instead of compressing and storing the C0 2> exhausting it as required directly through 
an Ohmjet which is designed to operate with an inlet pressure of one atmosphere. 

Systems were sized for both these methods and the component weights, peak power, 
and average power requirements were estimated for two recovery duty cycles. 

For all design studies in this section it is assumed that the space station has two 
identical Ohmjet systems, each operating on a 90 minute cycle with the cycles staggered 
in time by 45 minutes. Therefore, one system is applying thrust impulses for 45 mmutes 

while the other system is absorbing C0 2 . 

B. Description of Carbon Dioxide Recovery Systems - 

(1) General: 

The schematics of the two methods of C0 2 recovery systems are shown in Figures 
HI - 1 an d m-2. Continuous operation of a pair of these systems, located at diametrically 
opposite positions on the station is required for proper attitude control and orbit mainte- 
nance. During the first half of an orbit, one system (system A) is adsorbing CC> 2 from the 
cabin air, while the second system (system B) is going through its regeneration and re- 
covery cycle. During the second half of the orbit, system B is adsorbing CC> 2 while the 
sieve bed of system A is being regenerated and C0 2 recovered. Compensation for aero- 
dynamic drag on the space station is achieved if the Ohmjet of System A is directed to 
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SCHEMATIC DIAGRAM OF CO RECOVERY SYSTEM FOR DIRECT USE OF 

CO FROM MOLECULAR SIEVE BED, METHOD II 
2 

FIGURE III - 2 


2 




apply thrust toward the sun while the station is on the half-orbit approaching the sun, and 
the Ohmjet of System B applies thrust away from the sun while the station moves away from 
the sun, as shown schematically in Figure III 3. 

(2) Molecular Sieve Bed: 

As can be seen from Figures 111-1 and III-2, the molecular sieve bed and its elec- 
trical heater for regeneration are identical for both methods. The size and operation of 

this system will be discussed in detail in this section. 

The adsorbent material used in the C0 2 removal bed is Linde Type 5A molecular 
sieve material. During the adsorption cycle, all valves of the system are closed except 
V and V and predried cabin air is passed through the pelletized bed material, the tem- 
perature If which is approximately 77° during this phase of the cycle. At this temperature 
the molecular sieve bed material (See Figure III-4) wiU hold 0.08 pounds of COj per pound 
of adsorbent at a CO., partial pressure of 8 mm of mercury (the maximum permissible in 
the cabin air). Since it is required to recover the CX> 2 adsorbed by the bed, the normal 
method of regenerating the molecular sieve material by venting to space vacuum cannot be 
used. Instead it is proposed to heat the bed electrically to 392° F. during which time all 
valves would be closed. At this temperature the partial pressure of the CC> 2 is increased 
to one atmosphere absolute. Carbon dioxide can then be drawn from the bed at 760 mm of 
mercury. (This gas would be piped to the compressor inlet in Method I and directly to 
the Ohmjet in Method II.) The molecular sieve bed material at this temperature (392 F.) 

and pressure (760 mm of Hg) will hold 0.023 lbs. of COj per pound at adsorbent. Each 
regeneration cycle could therefore remove 0.08-0.023 = 0.057 lbs. of 0O 2 per pound of 

adsorbent . 

The energy required to heat each pound of the bed material for regeneration is 
CJT where C is the specific heat of the bed material and AT is the temperature rise. If 
we assume c/. 0.2 BTUAb. -°F and set AT = (392° F. - 77° F.) = 315°F. , then the energy 
required to heat the bed for regeneration is Q = 63 BTUAb. or 18.4 watt hoursAb. 

Upon completion of regeneration it is necessary to cool the molecular sieve bed 
down to its original temperature of 77° F. before the adsorption can begin. This means 
that the 63 BTU per pound of bed material must be removed. For this cooling cycle all 
valves except V3 are closed. The circulatory fan is started, opening check valve W & . 

The fan circulates the gas in this closed system through the bed to a radiant cooler and 
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thence back through the bed. Heat is removed from the circulating gas by the cooler and 
radiated to space. As the bed cools, the pressure in this closed loop drops. To prevent 
this pressure from dropping too low and decreasing the heat transfer capacity of the sys- 
tem, a small quantity of cabin air is sucked into the system through Valve V g and check 

Valve V to maintain the pressure near cabin ambient. When tae bed temperature is again 
4 

reduced to 77° F. the bed is ready to start another adsorption cycle. 

(3) Method I: Compressor, Storage Tank and Ohmjet: 

' ' o 

When the molecular sieve bed reaches the regeneration temperature of 392 F. 

Valve V 6 is opened and the centrifugal compressor is started Carbon dioxide is driven 

from the bed and enters the compressor at one atmosphere ar.d 392 F. It is compressed 

and is discharged from the compressor at 10 atmospheres and about 1000 F. through 

check Valve V^ into the storage tank. In order to purge the bed properly, a small quantity 

of C0 2 is bled through Valve V g and recirculated through the bed. The storage tank has 

the capacity to store the amount of CC> 2 recovered during one cycle. This storage tank 

acts as an accumulator which assures that C0 2 at 6 to 7 atmospheres is available to the 

Ohmjet which is operated for 2 seconds out of every 20 seconds during the regeneration 

cycle of 4.5 to 8.0 minutes each half -orbit. The complete system including the molecular 

sieve bed must be well insulated to minimize heat losses. Because the inlet pressure is 

between 6 and 7 atmospheres the size of the Ohmjet for Method I is relatively small. Also 

the heating element need only raise the gas temperature from 1000° F. to 1600° F. , which 

is the desired nozzle inlet temperature . 

(4) Method II: Direct Ohmjet, Operating at One Atmosphere Pressure: 

In this method, the compressor and storage tank are eliminated and the gas is 
supplied at one atmosphere and 392° F. directly to the Ohmjet through Valve V fi . In this 
method, there is no need for an accumulator, because the gas volume in the molecular 
sieve bed and piping is sufficient to maintain relatively constant pressure to the Ohmjet 
during its two seconds of operation every twenty seconds for 4.5 to 8.0 minutes each half- 
orbit. This complete system is also well insulated to minimize heat losses. With this 
method, the inlet pressure is low (one atmosphere), the Ohmjet is relatively larger and 
the heating element must be capable of raising the gas temperature from 392 to 1600 F. 
The inlet valve to the Ohmjet, however, operates at only 392° F. instead of at 1000 F. as 
int Method I and should therefore be more reliable. 
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C. Sizing of System Components - 

(1) Basic Assumptions: The sizing of each component for both methods is based on 
the following assumptions: 

(a) There are 21 men on the space station. 

(b) Each man generates an average of 1.88 lb. of CO- per day. 

(c) A maximum of 300 lb. -sec. impulse per orbit is required from the 
Ohmjets to control attitude and maintain orbit altitude of the station. 

(d) Station rotates at 3 rpm. 

(e) Ohm jet electrical efficiency is 90 per cent. 

(f) Nozzle thrust coefficient of Ohmjet is 1.7. 

(g) Specific impulse of CO^ is 60 sec. at 520° R. and varies as the square 
root of nozzle inlet temperature . 

(h) Cooling motor -fan is 50 per cent efficient. 

(i) Compressor has an efficiency of 70 per cent. 

(j) Compressor motor efficiency is 75 per cent. 

(k) Each kilowatt of power (average) requires 350 lbs. of auxiliary power 
equipment (e.g. solar cells and batteries). 

(2) Duty Cycles: 

The size and power requirements of each method depends directly on the duty cycle 
of the recovery system. In order to show the general effect of duty cycle on these basic 
parameters, the size, weight and power requirements of the components are derived for 
the following two cycles, based on a 90 minute orbital period (only cycle 2 has a longer bed 
heat up time and regeneration time to reduce peak electrical power demand): 


System A 

System B 

Duty Cycle 1 

Duty Cycle 2 

Adsorb 

L 

Heat Bed 

0 to 4.5 min. 

0 to 8 min . 



Regenerate 
and Pulse 

4.5 to 9.0 min. 

8 to 16 min. 


r 

Cool Bed 

9.0 to 45.0 min . 

16 to 45 min. 

Heat Bed 

Adsorb 

45 to 49.5 min. 

45 to 53 min. 

Regenerate 
and Pulse 



49.5 to 54.0 min. 

53 to 61 min. 

Cool Bed 

> 

f 

54 to 90 min. 

61 to 90 min. 
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(3) Molecular Sieve and CC^ Recovery System: 

(a) Adsorbent Required. The CC^ bed of each system must absorb all the 
CO„ generated by the 21 man crew during half an orbit (45 minutes). If we assume the 
average rate of CO 2 generation per man is 1.88 lb. /day, the total CO 2 to be adsorbed is 
39 lb. /day. During the adsorption cycle, each bed must therefore remove 1.24 lb. per 
cycle. As calculated previously, each pound of adsorbent can remove 0.057 lb. of CO 2 
per cycle. The required weight of adsorbent per bed is then 21.8 pounds. 

(b) Electrical Heating Requirements. The energy required to heat each 
pound of adsorbent from 77° to 392° F. was previously calculated to be 18.4 watt hours. 
The total energy to heat the required 21.8 lb. of adsorbent material is then 400 watt hours. 
This energy must be supplied during the bed heating cycle. The power level required to 
heat the bed is then 3000 watts for Duty Cycle 2. In addition it is assumed that 400 watts 
is required to maintain the bed at temperature during the regeneration cycle. 

(c) Cooling system . If we neglect the heat loss from the bed through the 
insulation, the amount of energy added to the bed during heating (400 watt-hrs.) must be 
removed from the bed during the cooling cycle to return the bed temperature from 392° F. 
to 77° F. The cooling system for the molecular sieve bed consists of a radiator and circu- 
lating fan in series with the molecular sieve bed. 

We define: w = circulation gas flow rate 
S 

T.^ = gas inlet temperature to radiator 

T o ut = g as outlet temperature from radiator 
T = average temperature of bed 

D 

T = average temperature of radiator 
R 

= radiation heat transfer rate 
A^ = radiator area 

If we assume T. = T„ 
in B 

T. + T 

^ in out 

T R= 2 

and Q r = «A r T r 4 

where a = Stefan Boltzman constant 

e = radiator surface emissivity ( 0.8) 
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where C p = 0.205 BTU/lb. - R for CC> 2 


C =0.20 BTU/lb. -°R for sieve bed 
B 


W = 21.8 lb. 
B 


From the above two equations, the molecular sieve temperature as a function of 
time is calculated and presented graphically in Figure HI-5 for two circulation rates and 

two radiator areas. 

The temperature histories show in Figure Ill-S do not account for any additional 
heat losses through the bed, piping or radiator internal wall, which would increase the rate 
of cooling. With this in mind, it is judged that a 10 square foot radiator w.th a 1. 17 lb. /mm. 
(~10 CFM) circulatory fan would be sufficient for the 36 minute cooling cycle of Duty Cycle 
1 and .ha, a 12 square foot radiator with a 1.75 lb. /min. (~15 CFM) circulating fan would 

be adequate for the 29 minute cooling cycle of Duty Cycle 2. 

The radiator is assumed to be a plate of 0.0625" aluminum with internal cooling 
channels. The inner surface is covered by one inch of insulation to minimize heat transfer 

into the cabin. 
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It is estimated that the radiator would weigh 1.69 lb. /ft. , consisting of 1.44 lb. 
of aluminum and 0.25 lb. of insulation. Duty Cycle 1 requires a 10 ft. of radiator which 
would weigh 16.9 lbs. and Duty Cycle 2 requires a 12 ft. radiator which would weigh 


20.3 lbs. 

The flow rates of the circulating fan have been established and an estimate of the 
head against which it must operate is required to establish its power requirements. Con- 
servative estimates of the head requirements are tabulated below for 10 CFM and 15 CFM 

flow rates: 


ITEM 

Velocity head for 2" i.d. pipe 
Check valve pressure loss 
Radiator pressure loss 
Pipe Friction 
Bed loss 

TOTAL HEAD 

Air Power 

Efficiency of fan 

Efficiency of motor 

Power requirement of motor 


10 CFM 

0.021 in of H 2 0 
0. 100 in of H 2 0 
1.400 in of H 2 0 
0.063 in of H 2 0 
1.400 in of H 2 0 

2.984 in H 2 0 
155 ft. lb. /min. 
50% 

70% 

443 ft. lb. /min. 
or 10 watts 


15 CFM 

0.046 in of H 2 0 
0. 100 in of H 2 0 
2 . 200 in of H 2 0 
0. 142 in of H 2 0 
3. 150 in of H 2 0 

5.638 in H 2 0 
4.38 ft. lb. /min. 
50% 

70% 

1250 ft. lb. /min. 
or 28.3 watts 


(d) Compressor - 

The compressor in Method I must pump the 1.24 lb. of CC> 2 generated from 
1 to 10 atmospheres during the regeneration cycle. The work required to compress this 
gas isentropically is expressed by 



where Ah = work per unit weight of gas 
s 

R = gas constant 
T = inlet temperature 
k = ratio of specific heats 
P = inlet pressure 
P = discharge pressure 


£0 


Substituting the values R = 35.1 ft. /°R., ^= 852^., P^Pj = 10, and k = 1 . 28 in the above 

expression we get Ah - 70, 000 ft. The energy required to isentropicaliy compress the C0 2 

s 

is therefore 86, 900 ft. lb. /cycle. Assuming adiabatic compression with a compressor of 

o 

efficiency of 70 per cent, temperature rise in the compressor would be 950 F. The com- 
pressor discharge temperature would be 1800 R. or 1340 F. Allowing for heat conduction 
into the compressor during the short period of compressor operation, it is probable that the 
gas exit temperature would be reduced to an actual value near 1000 F. 

The power to drive the compressor should be supplied at a high speed (24, 000 rpm 
motor). If this motor has 75 per cent efficiency, the input power required by this 
motor is 840 watts for Duty Cycle 1 and 470 watts for Duty Cycle 2. Motors of this size 
operating at 24, 000 RPM weigh about 2-1/2 to 3 lbs. per kilowatt. Therefore, the motor 
weight is estimated at 2.2 lbs. for Duty Cycle 1 and 1.4 lbs. for Duty Cycle 2. 

A high speed 2 or 3 stage compressor for such an application would have a diameter 
of about 10" and width of 1.5" with a metal solidity of about 1/3. Its weight therefore is 
estimated as: 

tt (lO)^ x j x jyg x q 286 = 11.2 lbs. 

T 

There would be little difference in weight in the compressor required for either Duty Cycle. 

( e ) Storage Tank . The storage tank is designed to hold the complete amount 
of CO regenerated in one cycle at a temperature of 1500° R. and a pressure of 147 psia. 

2 3 

The tank volume is calculated to be 3 ft. . The lightest structure which will hold this volume 
is a spherical tank with an internal diameter of 1.79 ft. The required wall thickness is 
calculated to be 0. 020 in. for a working stress of 40, 000 psi in stainless steel. The pressure 
shell of the tank would weigh 8 . 3 lbs . 

To minimize heat loss from the stored gas, the tank is covered by 2" of insulation 
covered with an external cover of 0.010" thick aluminum. The weights of these two items 
are 6 lbs. and 2 lbs. respectively making the total weight of the storage tank 16.3 lbs. 

(f) Ohmjets. During each orbit, 2.48 lbs. of CC> 2 are exhausted through the 
Ohm jets . In order to obtain the total 300 lb-sec/orbit impulse required, the specific im- 
pulse of the CO must be 121 seconds. For CO this value of specific impulse corresponds 
to a gas temperature of approximately 1600 F. 

During the regeneration cycle, the Ohmjet is given a 2 second pulse once every 20 
seconds (once each revolution of the station). The heating element of the Ohmjet is operated 
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for 1-1/2 seconds preceding and during the 2 second pulse of gas flow for a total of 3-1/2 
seconds each 20 seconds. 

For Method I, the total energy put into the gas by the Ohmjet is 90 watt hours/orbit. 

If we assume an Ohmjet efficiency of 90 per cent, the energy supplied to the Ohmjet is 100 
watt hours /orbit. 

For Duty Cycle 1, there are 27 pulses per orbit, therefore each pulse must consume 
3.71 watt hours in 3-1/2 seconds. The power required by the Ohmjet is therefore 3810 watts. 
For Duty Cycle 2, there are 48 pulses per orbit, thus the power required by the Ohmjet is 
2140 watts. 

A total impulse of 300 lb. -sec. /orbit is assumed to be required rom the Ohmjets. 

This impulse equals the product of the thrust, F, and the time of operation per orbit. The 
Ohmjet thrust is also expressed by 

F ■ Wp 

where P = the Ohmjet nozzle inlet pressure 

= the Ohmjet nozzle throat area 

C = the nozzle thrust coefficient 
F 

Assuming P = 100 psia, and C = 1.7: 

VA square inches 

For Duty Cycle 1, the Ohmjet thrust is 5.56 lbs. , which requires a nozzle throat 
2 

area of 0.0328 in (0.204 in diameter). 

For Duty Cycle 2, the Ohmjet thrust is 3.13 lbs. , which requires a nozzle throat 
2 

area of 0.0184 in (0. 153 in diameter). 

Using the same procedure as above, the Ohmjet energy requirements for Method II 

are calculated to be 201 watt hours /orbit. The power level required by the Ohmjet is 7690 

watts for Duty Cycle 1 and 4310 watts for Duty Cycle 2. For Duty Cycle 1, the required 

2 

Ohmjet thrust of 5.56 lbs. requires a nozzle throat area of 0.223 in (0.532 in diameter) 

and for Duty Cycle 2, the Ohmjet thrust of 3. 13 lbs. requires a nozzle throat area of 
2 

0. 125 in (0.400 in diameter). 

D. Electrical Power Requirements - 

The power requirements for each method and both duty cycles are summarized in 
Table I and shown graphically as a function of time in Figures III-6, 7, 8, and 9. 
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TABLE HI - 1 


ELECTRICAL POWER REQUIREMENTS FOR OHMJET FOR TWO DUTY CYCLES 
DUTY CYCLE 1 


- TIME 

FUNCTION 

METH 
Peak Power 
(watts) 

[ODI 

Total Energy 
(watt-hours) 

METHOD 
Peak Power 
(watts) 


0 *o 4 .5 min 

Heat bed B 

5330 

400 

5330 

400 

4 .5 to 9 min 

Regenerate and Pulse 

5050 

143 

8090 

131 

9 *"0 45 

Cool bed B 

10 

6 

10 

6 

4„ to 49-1/2 

Heat bed A 

5330 

400 

5330 

400 

4° -1/2 to 54 

Regenerate and Pulse 

5050 

143 

8090 

131 

5^ to 90 

Cool bed A 

10 

6 

10 

6 

1 rbit 

TOTAL 

5300 

Average Pow 

1098 

:r = 732 watts 

8090 

Average Powej 

1074 

= 716 watts 


DUTY CYCLE 2 


"TIME 

FUNCTION 

METt 
Peak Power 
(watts) 

IOD I 

Total Energy 
(watt -hours) 

METHOD 
Peak Power 
(watts) 

>11 

Total Energy 
(watt-hours) 

1 t 8 min 

Heat bed B 

3000 

400 

3000 


to 16 min 

Regenerate and Pulse 

3012 

166 

4710 1 

154 

6 5 48 min 

Cool bed B 

28 

13-1/2 

28 

13-1/2 

5 to 53 min 

Heat bed A 

2000 

400 

3000 


3 > 61 min 

Regenerate and Pulse 

3012 

166 

4710 

166 

1 to 90 min 

Cool bed A 

28 

13-1/2 

28 

13-1/2 

i 

o 1 

TOTAL 

3012 

Average Power 

1159 

= 772 watts 

4710 

Average Powei 

1135 

= 756 watts 
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Although the average power requirement of the system is approximately 750 watts, 
higher peak power is required to heat the molecular sieve bed and pulse the Ohmjets . It 
should be noted, however, that all these peaks occur when the station is in the sun and 
therefore can come directly from the solar cells and not from the batteries. Since these 
peak loads do not exceed the normal excess power generated by the solar cells for charg- 
ing the batteries, the maximum demands of the Ohmjet can be supplied simply by a tempo- 
rary reduction of the charging rate of the batteries. 

E. Weight Estimates - 

The weight estimates for the components required for complete Ohmjet systems are 
given in Table IIT-2. 

The total system weights include all components including additional electrical 
power system weight needed for the operation of one complete space station system, and 
does not include any extra components or subsystems which would be required for redun- 
dancy. Even if another complete system were added for 100 per cent redundancy, an Ohmjet 
system shows considerable weight savings compared to the 5840 lbs. /year of bipropellant 
fuel (I = 300 seconds) which would be required to do the same job. 

Although Method II offers considerable weight savings as compared to Method I, 
for both Duty Cycles, the weight difference is not so decisive as to rule out Method I. 
However, factors such as reliability and simplicity make the use of Method II worthy of 
more detailed investigation. 

F. Discussion - 

Both methods of utilizing the waste COj of the space station for attitude control In- 
vestigated in this report are feasible and workable systems . Each method has certain 
advantages and disadvantages. These are listed below: 

Method I - Compressing and storing C0 2 prior to use in the Ohmjet 

Advantages 

(1) Smaller Ohmjet thrust unit required because of higher inlet pressure. 

(2) Smaller temperature rise required in the Ohmjet because of high compressor 
discharge temperature. 

(3) Storage tank acts as an accumulator to smooth over transient peaks anjd varia- 
tions in C0 2 generation rate. 


TABLE 111-2, WEIGHT ESTIMATES FOR OHMJET SYSTEMS 

Weights in Rounds 

— “ METHOD I 1 METHOD II 

ITEM putyCv cle T I Du ty Cvcl eljDutv Cycle_L 

Molecular Sieve Bed N r * N , c . | N.C. f 

Bed Tank — 1 fl N " jv 

Bed Material (21 .8# Linde Type 5A) - 1 o.*50£ 0*40# I 0.50# C 

Heater - 1 2 56 1 2.56 2 

Insulation ( 1'* of 3#/cu ft Foam) V~T90 1 

Insulation cover ( .010" Aluminum) 1 — : L I - 


N.C. 

N. C. 

O. 50* 
2.56 
1.90 


N.C. 

N. C. 

O. 40# 
2.56 
1.90 


TOTAL Molecular Sieve Bed 


nnlpr Plate 10 .080" Aluminum 


4.96 


4.86 


4.96 


4.86 


insulation Cover ( .010" Aluminum 


TOTAL Radiator 


Circulating System 


Circulating Motor- Fan 


V9 - 0.1” of HoO check valve 


and Insulation 



TOTAL, Valves and Remaining Piping 


* N .C . - Denotes not chargeable to Ohmjet because 

item is required for life support . 3“ 











































TABLE IU-2, WEIGHT ESTIMATES FOR OHMJET SYSTEMS (CONT*D) 














Disadvantages 


(1) Heavier. 

(2) Complex, requires more components such as compressor, storage tank, 
check valve and pressure regulator. 

(3) Ohmjet control valves must handle high pressure (150 psia) and high temperature 
(1000° F.) gas. 

Method II - Direct use of regenerated CC> 2 from molecular sieve bed into the Ohmjet. 
Advantages 

(1) Simple, requires only a few basic components. 

(2) Lighter weight. 

(3) Ohmjet valve must handle only low pressure (14.7 psia) and low temperature 
(392° F.) gas. 

(4) The functual simplicity of Method II should lead to a higher reliability than that 
of Method I. 

It is concluded that:: 

(1) The Ohmjet thrust system is q simple, feasible and practical system for obtain- 
ing attitude control by the use of waste CQ 2 from a manned space station. 

(2) Direct use of C0 2 during regeneration of a molecular sieve bed provides a 
simplex* lighter, and potentially more reliable system than one using a CO, compressor 
and storage tank. 
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IV design studies of a solar-heated pulse 

ROCKET (HEPOJET) THRUST SYSTEM 


A . Design Requirements 

•me design concept of the Heliojet is to utilize a paraboloidal solar collector 
with a cavity receiver containing a l)eat exchanger, consisting of helical coils of 
metallic tubing. The unit is mounted on the sunward side of the space station so 
that the solar energy collected by the mirror is stored as thermal internal energy 
in the metallic heat exchanger tubes . Periodically, die propellant gas (in this case 

C0 2 ) is passed through the heat exchanger tubes where it is heated and then expanded 
through a nozzle to produce a thrust impulse. 

A basic advantage of the Hellojet, as compared to the Ohmjet. stems from 
potentially higher efficiency due to the fact that solar energy is used to raise the 
propellant temperature directly, rather than through die intermediate process of 
rting solar energy to electrical energy . For the design study of the Hellojet , 
the same inlet conditions of the gas were assumed as for the Ohmjet and the objective 

has been to determine whether a realistic design configuration can be obtained with 
specific impulse comparable to that of the Ohmjet . 

Based upon the stud.es in Section 111, d,e following conditions were assumed for 
the C0 2 at the inlet to the Heliojet heat exchanger tubes: 


Temperature 100°F 

Pressure 150psia 

Mach Number 0.05 


The above gives a gas velocity of 45 ft/sec and a density of 1 .1 lbs/ft 3 . The 
corresponding mass velocity, G, is 49.5 lbs/ft 2 -sec. 

It was established that if the <X> 2 production from 21 men were heated to 1<S00°F 
the average rate a, which impulse could be applied to die station would correspond te ' 
a continuous dirust of 0 .0547 lbs . This produces two additional Heliojet design specifi- 
cation.: (1) me average gas temperature during the impulse must be about 1600°F and 
2, me average gas flow rate (Ib/sec) must be such diet an average thrust value of 0.0547 lbs 
■s produced. It is to be noted that die latter specification is based upon die bluest 
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rate of thrust application , i .e . when all of the available CC> 2 is used for this purpose , 
which will occur during periods of maximum gravity gradient torque . 

It is clear that the impulse duration is of primary importance in determining the 
average gas temperature attainable . Also , allowance must be made for the occurrence 
of eclipse during periods of maximum required thrust . Thus if we assume maximum 
eclipse time to be 35 minutes for an orbital period of 95 minutes , the thrust level 
must be increased by a factor of 1.585 so as to make up for the time in the dark when 

the Heliojet cannot operate . 

Since the attitude control pulse must be applied at a predetermined angular 
orientation, the frequency of thrust application is related to the rate with which the 
station rotates about the sun-pointing axis . It will be assumed that this spin rate 
is 3 RPM . Hence , an impulse could be applied once every 20 seconds or at multiples 

of 20 second intervals , i .e . , once every 40 sec . , 60 sec . , etc . 

If thrust were to be applied once per rotation , the required impulse would 
Ft = (0 .0547) (20) ( 1 .585) = 1 .735 lb-sec per rotation ( 1) 

It will be assumed that the design thrust level is 1 .735 lbs . This permits a 
choise of t , the thrust duration , depending upon the frequency of thrust application . 
Thus , if the impulse is applied once every rotation, t=l sec . , if applied every 40 
seconds , t = 2 sec . A thrust duration much longer than 2 seconds is not as efficient 
in that the average moment arm of the thrust is reduced as the angular firing angle 
exceeds 36° . 

The various requirements which the Heliojet design must meet can be further 
defined as follows: 

(a) Since allowance should be made for prolonged periods of time in sunlight 
without thrust application , the maximum heat exchanger temperature must be 
kept to a reasonable value so as not to damage the unit . If we assume the tubes 
and cavity liner to be made of Niobium , thus temperature should not exceed 
2500°F . This requirement is used to limit the concentration ratio of the solar 

collector . 

(b) During C0 2 flow, the gas temperature at the exit of the heat exchanger 
tubes must average 1600°F . in order to produce the desired specific impulse of 
120 sec . This requirement is the principal consideration in the analysis of 
transient heat-transfer from the tube walls to the gas . 
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(c) In the time available between thrust application the tube wall temperature 
must rise to an intial temperature high enough to be consistent with (b) above . 

(d) Collector and cavity geometry must permit angular misalignments with 
respect to the sun of one degree or more without affecting operation of the 
system . 

In order to define the design of a Heliojet thrust system , analyses were 
performed in three separate areas: optical design, receiver heat-up rates , and 
transient heat-transfer to the gas . These are described below, followed by a 
description of a configuration meeting the various requirements defined . 


B . Optical Design Considerations 
( 1) Concentration Ratio 

Since the exterior of the cavity will be insulated to minimize heat loss 
from the walls , and since precautions will be taken in the design to limit losses 
due to direct reflection from the cavity opening to less than 2 per cent , the heat 


’-’dance equation can be written as: 

QsAp-AjTT^WCpf: 

where 

Q s 


( 2 ) 


T 

cr 


Solar heat flux = 442 BTU/hr-fr 

Area of paraboloid 

Area of cavity opening 

Heat-exchanger tube wall temperature , C R 

Stefan Boltzman constant = 0.174 x 10~ 8 

weight of receiver tube walls , lbs . 


BTU 

hr -ft'* (°R) 4 


i .e . 


W 

Cp - specific heat of tube material, BTU/lb. -°F . 

The maximum receiver temperature will be reached in the steady-state , 
when dT/dt = 0 . Hence 

T =f <9aK*B> 1 1/4 <3) 

max L(<T)(A c )J 

Since A p /A c defines the concentration ratio, C, its value can be determined 
for a given value of T 

max • 

^ or T max - ^ . (2,960 R), the maximum allowable concentration ratio 

is: 


p = C = 300 


(4) 
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(2) Losses due to reflection - 

In order to approximate a black-body cavity absorber , the amount of 
radiation passing directly out of the cavity opening after reflection from interior 
surfaces must be minimized . This is accomplished by making the ratio of cavity 
aperture area A to cavity internal surface area A. as small as possible.. As 
shown in Reference IV- 1, when this ratio is of the order of 0.02, the fraction 
of light passing directly out of the opening would be about 2 percent. This 
fraction is , for all practical purposes , independent of cavity geometry , i .e . 
whether conical , spherical or cylindrical in shape . It is also not appreciably 
affected by the value of interior surface emittance if the latter is in the range 
of 0 .3 to 1.0 over all wave lengths . 

In the Heliojet design , the total internal area of the cavity is made up of 
the cavity wall and the exterior of the heat exchanger tubes . It is assumed that 
these tubes will be spaced apart by one tube diameter in order to permit reflected 
radiation within the cavity to heat the tube wall uniformly on all sides . 



Figure IV- 1 Geometry of Heat Exchanger Tube 
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Since the cavity diameter is approximate! equal to the coil diameter d^ 
A i =n2Dd h N + 2nDd h N (5) 

where N = number of helical turns 
The first term above represents the tube heat transfer area A 


Hence 


HT 


A i = A ht + TT A ut = 1 -64 A IirT , 
i MI tt HT HT 


In order that 


( 6 ) 


c < 

— i 0.02 

1 

It is necessary that 
A 


_c_ (1.64)(0.02) = 0.033 


HT 


For design purposes a value of A c /A HT =0.03 will be used . This will 
maintain losses due to direct reflection to less than 2 percent (Ref . IV- 1). 

(3) Angular misalignments - 

Reference 1V-2 provides the basic relationships and criteria for estimating 

the amount of angular misalignment which can be tolerated for a given sire cavity 

For nomenclature as shown in Figure 2, pertinent relationships are summarized 
in Fig. 1V-2 . 


P 

f 

P 


r 

a 


2f 

r+ees 4L 


= D p < 1+C0S © R ) 


4 sin ©, 


2 sin © 

K. 

0.00931 f 


( 1 + cos©^ COS€^ 


s = " f a Dp (3-cosefr) 

*T COS ^R 2 sin 20 R 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 
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The distance S in equation (12) is the maximum motion, relative to the 
center of the circular solar image , of a point in the focal plane reflecting from 
the rim of the mirror due to a misalignment angle o< with respect to the sun . 

To minimize motion of the image it is desirable to choose a value for the rim 
angle . © r > at which 

d(S><) _ 

The corresponding value of the rim angle is ©^ = 42° . This value is 
therefore used to define the geometry of the paraboloid. 

According to Reference IV-2, the energy distribution within the aberrated 
solar image can' be approximated by assuming it to be of constant intensity within 
the circular solar image and then falling off linearly until it reaches zero intensity 
at r = r^ . On this basis , the image which is to be confined within the cavity 

aperture has a radius r & given by equation (11). In terms of the diameter EL and 

o * 

for © = 42 , this becomes - 
r 

r fl = 0.00467 $13) 

The angle o< at which r & will reach the edge of the cavity aperture and 
still be confined within it is found from 


r = S + r 
c a 

(14) 

which for ©^ = 42° yields 


C* =50 .6/ i£.\- 0 * 236 (degrees) 

(15) 


Having fixed the concentration ratio at 300, the ratio 
also becomes fixed . Thus 



in equation ( 15) 


r c = _L 
Dp 34.7 


and the misalignment angle which can be accepted without affecting operation of 
the solar heating system is 

0< =1.22 degrees 

max ® 

C . Receiver Heat-Up Rates - 

The rate at which the heat-exchanger wall temperature will rise during the 
time when there is no propellant flow can be calculated from equation (2): 
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dT A p , _ (£ /Q S ) T 4 

dt = p 6 C p X ht (A p /A c ) (18) 

The concentration ratio i^/A c has been previously selected to be 300. The ratio 
A^/Atj T has been selected to be 0.03. Hence 

C n X 

A 300 A 

P- = — — = (0.03) (300) = 9 (19) 

a ht a ht 

For tubes made out of Niobium to withstand the maximum temperature of 2, 500 F, 

C =0.1 BTU/lb — °F 
P 

o = 524 lbs /ft 3 
Equation (18) becomes 

r—^ 4 ] <*» 

m U* — 

where 6 = wall thickness in thousands of an inch (mils) 

m 

For the present purpose, an explicit integration of equation (20) is not warranted. 

To obtain the variation of wall temperature with time for assumed initial temperatures T 

dT* ^ 

and wall thickness 6 , a step-wise numerical process was used. It was assumed-nj- is 

' m . m 

constant over temperature intervals of + 50 F. around the value of T for which dtT is 
computed. Results are plotted in Figure IV-3 for assumed initial temperatures o?750° F., 
950° F. and 1, 250° F. For each of these values, the temperature rise is shown for wall 
thickness of 5 and 10 mils. 

It is evident from Figure IV-3 that if the wall temperature is to rise to a value 
substantially above 1,600° F. in one or two revolutions of the station, the initial tempera- 
ture (at the end of a gas flow pulse) must be in the order of 1,250° F. If operation of the 
Heliojet is to occur once every rotation of the station, i.e. every 20 seconds, a wall 
thickness of frlfejils ’Adll.fee required to reach a temperature of about 1, 950° F. Similarly, 
if thrust is to be applied every two rotations, the same temperature could be reached with 
a wall thickness of 10 mils. 

Without increasing the thrust level beyond the assumed value of 1.735 lbs. , the 
greatest flexibility in programming the frequency of applied impulses will be obtained if 
it is possible to actuals the Hello jet every rotation of the station. For thii reason a wall 
thickness of 5 mils is chosen. 
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D. Transient Heat-Transfer Analysis 

Appendix A contains a detailed description of the transient heat-transfer 
analysis and provides relationships for the variation of gas and wall temperature 
with time and as a function of the various design parameters . Figures A- 2 and 
A- 3 of the above appendix give dimensionless design curves for the gas and wall 
temperatures , respectively . These curves provide a means for selecting the 
diameter and length of the heat exchanger tubing and the thrust duration . The 
procedure is essentially one of trial and error, for the values ultimately chosen 
must meet the requirements for wall temperature described in the preceding section 
and must also match the design curves for both wall and gas temperature during 

transient heat-transfer . One set of values which meets these conditions has been 
found to be: 

Tube inside diameter - D = 0.10 inches 
Length of one tube - L = 30 inches 

Pulse duration - t = one second 

Referring to Figures A-2 and A-3 of Appendix A, the above design assumptions 

yield the following values for the coefficients K &K 5 

2 3 


K 2 = 4 S \%’ - 4 •« 

K = St (£)_ = 2 06 

3 <Cp w ) <p6) 2 ' 06 


( 21 ) 

( 22 ) 


The Stanton Number <St> has been evaluated in accordance with the relationships 

defined in Appendix A; all other parameters have been selected in the preceding sections . 
Using the above: 


T wa<*> ’ 100 


= 0.63 


(23) 


of 


wa 1970 - 100 

This yields an average wall temperature at the end of the one second impulse 


T wa (1) = 1 ’ 275 F * 

This value, and the initial value of T^ = 1970°F . used in equation < 23) , match the 

initial and final conditions of temperature rise during the tube-wall cycle (See Figure 
IVtJ for 6 = 5 mils). 
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The highest gas exit temperature , reached at the start of me impulse , is found 


from the value of at t = 0 . This yields 


Hence , 


-£T~ = T (0) - 100 = 1 “ e - 0 * 988 

ge ge 

1970 - 100 


T (0) = 1945 F . 

ge 


(24) 


Exit gas temperature at the end of the one second impulse is found from 


T (1)-100 


*y i)L r 


100 


= 0.625 


(25) 


T (1) = 1270 F 
ge 


(It is to be noted that this value of gas temperature cannot be compared directly 
with the average wall temperature of 1,275° F . since the latter represents an average 

over the tube length . See Appendix A for further discussion .) 

The average value of gas exit temperature during the impulse is therefore 

1945 + 1270 = 1607° F 

2 

The average value of gas exit temperature needed to produce the assumed specific 
impulse of 120 seconds is 1,600°F . 

E . Heliojet Configuration 

The preceding analysis Establishes the length and diameter of the tubing m the 
Heliojet heat exchanger: L/D = 300 , D = 0 . 1" , L = 30” . However . the diameter thus 
chosen will not accomodate the required mass flow in one tube . The number of parallel 

tubes needed is found from 

a f =nJd 2 =!l (26) 

F * G 

where N = number of tubes 

.incew - j£ = 4ir = 0.01448 U>/sec < 27 > 

4 8p 

and G = 49 .5 lb/ft 2 - sec 
The above yields 


N = 5.36 



To produce a round number of tubes, say N = 5, the diameter should be 0.1035" 
rather than the value of 0 . 100" previously used . This small change does not appreciably 
affect the preceding calculations because D enters through the Stanton number both into 
the calculations of K 2 as well as . However , L/D must be maintained at 300 since 
it enters only into . To do this the length of each tube is increased to 

L = (300)(0 .1035) = 31.1 inches 

Using five tubes in parallel wound on a 3" diameter helix results in a heat 
exchanger as shown in Figure IV-4 . Also shown are related sizes of the concentrator 
and cavity based upon the design criteria described in preceding sections . 

F . Weight Estimates 

The estimated weight breakdown of the Heliojet thrust unit shown in Figure IV-3 
is shown in Table I . 


Table I - Estimated Weight Breakdown of Heliojet Thrust Unit 


Concentrator mirror (2 ft . diam .) 

Cavity Structure - 0.020" Niobium walls and radiation 
shielding 

Heat exchanger tubing 
Nozzle 

Concentrator support structure 

Total Thrust Unit Weight 


1 . 5 lbs . 

0 . 5 lbs . 
0.1 lbs. 
0 . 2 lbs . 
0.7 lbs. 

3.0 lbs. 


The weight of the CC> 2 recovery and storage equipment for a Heliojet thrust 
system is estimated on the basis of calculations made (see Section III) for an Ohmjet 
system , using Method I , i .e . using a compressor and storage tank to recover CO at 
a pressure of 10 atmospheres . Since the Heliojet uses C0 2 every 20 or 40 seconds 

,e the space station is in sunlight, it is necessary to store C0 2 during periods between 
regeneration of the molecular sieve beds . It is assumed that two sieve beds are used 
and each is regenerated once per orbit. It is conservatively assumed that heat losses 
will occur from the storage tanks so that the Heliojet inlet temperature is only 100°F . 
even though the compressor discharges CC> 2 into the storage tank at 1000°F . 

The estimated weight of components for a complete Heliojet system is shown in 
Table IV-2 . 
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Concentration Ratio : 300 
Internal Cavity Design: 

5 Tubes in parallel . Tube Material: Niobium . 

Tube ID: 0 . 104 inch . Wall Thickness: 0 .005 inch . 

Tubes are wound on a 3 .0 inch diameter helix with 
0.114 inch spacing between adjacent tubes . Developed 

length of each tube is 31 inches . Number of turns: 3 .4 per tube . Total no . of turns: 17 
CONFIGURATION OF HELIOJET THRUST UNIT 


FIGURE IV- 4 
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Table 1V-2 


Estimated Weight Breakdown of Heliojet Thrust System 


Two Solar-heated thrust units 

Two Heaters and Insulation for molecular 
sieve beds 

Two Radiators to cool sieve beds 

Two circulating systems to cool sieve beds 

Two compressor systems 

Two CO„ storage tanks 

Valves and piping 

Electrical controls for valves 

Total Heliojet System (no redundancy) 


Weight (lbs) 
6.0 


10.0 

40.0 
5.5 

28 .0 
32.6 

9.7 

6.0 

137 .8 lbs 


G. Conclusions 

1 . In order to compensate for the maximum anticipated gravity gradient 
torque by the use of a solar heated pulse rocket (Heliojet), it is desirable to design 
for a single thrust impulse to be delivered every one or two rotations of a space 
station while it is in sunlight. If corrections are made less frequently, the size of the 
solar collector would become excessive . The Heliojet system requires ttiat carbon 
dioxide be recovered , compressed and stored so that it may be used every 20 or 40 

seconds during periods of maximum torque requirement. 

2 . The Heliojet system analyzed above when operated on every rotation of the 
station produces the same performance as an Ohmjet system which is electrically 
heated to 1600°F . During periods of less (ban maximum torque requirement, the 
same Heliojet design may have a longer heat up time and would achieve higher average 
nozzle inlet temperatures . It has been estimated that if thrust is applied every two 
rotations , the average nozzle inlet temperature would be 1900°F . , leading to a 10 
percent increase in specific impulse to 130 seconds , as compared to 120 seconds 

for the reference design . 

3 . Operation of a Heliojet system requires that the space station spin vector be 

oriented to the sun with an accuracy of + 1.25 . 

4 . The estimated Heliojet unit weight of three pounds is low enough that two or more 

»«4r«i could be used for redundancy . 


55 



5 The estimated concentrator diameter of two feet is small enough to permit 
a simple one-piece mirror construction and die required concentration ratio of 300 
Will be ea- Pv achievable . (Concentration ratios of 1000 or more are required for 
solar-thermionic power systems) . 
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V EXPERIMENTAL TESTING OF OHM JET THRUST UNIT 


A. Test Objectives - 

The objectives of the experimental test program were to determine the thrust, 
specific impulse, and total impulse per pulse for an Ohmjet thrust unit operating with 
carbon dioxide as the propellant. During the test program, it became apparent that the 
nozzle inlet temperature was a sensitive function of the heat transfer coefficient from 
the helical heater wire to the gas, for which no prior empirical data were available . 
Therefore, preliminary measurements were made to determine the heat transfer coef- 
ficients, as well as the electrical efficiency for steady state operating conditions. 

B. Description of Ohmjet Unit - 

Figure V - 1 shows a drawing of the Ohmjet unit used in the test program . It 
consists of a tubular stainless steel chamber with a replaceable converging orifice as the 
nozzle, an internal helix of Nichrome wire as the heater element which is electrically and 
thermally insulated from the tube wall by a sheet of flexible asbestos paper, and electrical 
bus -bar connections to permit application of up to 220 volt, 60 cycle power to the heater 
element. 

The unit was designed for a OOj inlet pressure of 150 to 200 psia and to operate 
with heater element temperatures up to 1850° F. and nozzle inlet temperatures up to 
1600° F. The nozzle plate used in the experiments had a converging orifice with a threat 
diameter of 0.0625 inch which exhausted to the atmosphere. 

Carbon dioxide was supplied to the Ohmjet from a pressure tank by a pressure 
regulator and the flow was controlled at the inlet to the Ohmjet be a fast -acting solenoid 
valve (Maratta Valve Company Number MV- 100) which had been modified to attain an 
opening time and closing time of approximately 10 milliseconds. 

The heater coil consisted of a helix of 20 mil diameter Nichrome V wire with an 
outside coil diameter of 0. 180 inch and approximately 300 turns between the electrical 
contacts. The resistance of the heater at 70° F. was measured to be 22.55 ohms. There- 
fore, upon rapid application of 220 volts, the initial current was 9.75 amperes and the 

power, 2150 watts. The mass of the heater coil between contacts was calculated 
to be approximately 6.6 grams. Assuming a specific heat for Nichrome of 0. 11 cal/gm C, 
the initial heating rate with 220 volts would be: 
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OHMJET THRUST UNIT 
Figure V - 1 



<u _ p_ _ 2150(0.2388) _ 780 ° c / gec _ 1400 ° F/sec. 

S' WC (6.6)(0. 11) 

P 

C. Discussion of Ohmjet Testing under Pulsed Operating Conditions - 

The measurement of transient thrust levels below one pound with rise times of 
the order of 10 milliseconds poses severe problems in achieving adequate sensitivity 
and avoiding ringing and overshoot in recording the thrust vs time curves. A reasonably 

satisfactory method was achieved in the following manner. 

As shown in Figure V - 2 and V - 3, the Ohmjet thrust unit with its control valve 
attached was suspended on a thin elastically bent strip of stainless steel, 0.010 thick 
by 0.500" wide by 2.0" long which had strain gages (Baldwin Lima Hamilton Number 
C-14) mounted on each side. One strain gage was connected in a bridge circuit having a 
DC applied voltage of 90 volts. The output of the bridge was amplified in a high gain DC 
amplifier (Tektronix Model D) and applied to the vertical deflection plates of a Tektronix 

Model 53 1A oscilloscope. 

Carbon dioxide was supplied at a regulated pressure up to 150 psig through a long 
flexible tube to the inlet port of the control valve. The control valve was actuated by a 
24 volt DC power supply with opening and closing times programmed by a mechanical 
switch programmer (Eagle Model A6). When the control valve was opened, the thrust 
developed by the Ohmjet reduced the weight of the Ohmjet unit applied to the strain gage 
element and the thrust time curve was recorded on the oscilloscope as shown in Figures 
V - 4 and V - 5 . 

The strain gage system was calibrated statically by applying a vertical upward 
force on the Ohmjet with a spring balance and observing the resultant oscilloscope deflec- 
tion. The calibration curve is shown in Figure V - 6. 

With the stainless steel suspension strip initially bent to the degree necessary to 
obtain adequate bending stress in the strain gage to measure thrust levels of less than one 
pound, it was found that the naturally frequency of the suspension system for vertical 
vibrations was approximately 10 cps. The lO.cpe ringing oscillations shown in Figure V# 4 
were induced by the rapid application of thrust with insufficient system damping. In these 
tests a small piston and cylinder dashpot filled with air was connected to the Ohmjet to in- 
troduced damping by viscous and sHdfegifrtettaa.. However, with this design it was not 
possible to achieve a large fraction of critical damping, even with water in the dashpot. 
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OHMJET THRUST TEST STAND 
Figure V - 2 



PHOTOGRAPH OF OHMJET THRUST TEST STAND 
FIGURE V - 3 


& V 


Nozzle Throat Diameter = 1/ 16 in . 
Average Thrust = 0.31 pound 
Duration = 1.2 seconds 



THRUST - TIME RECORD AT 125 psig REGULATED PRESSURE 
FIGURE V - 4 



Nozzle Throat Diameter = 1/16 in. 
Average Thrust = 0.375 pound 
Duration = 1.2 seconds 


THRUST - TIME RECORD AT 150 psig REGULATED PRESSURE 


FIGURE V - 5 











Therefore, all oscillograms obtained exhibited some degree of initial overshoot and ringing. 
However, the average amplitude of the thrust and the total pulse duration were quite accu- 
rately repeatable from one pulse to the next. 

By comparing measured thrust levels with the thrust calculated on the basis of 
throat area, nozzle inlet pressure, and orifice thrust coefficient, it was estimated that the 
pressure drop from the regulator through the control valve and the Ohmjet heat exchanger 
section was approximately 28 per cent of the upstream regulated pressure for the test 
conditions of Figures V-4 and V-5. 

Transient tests were made with electrical power up to 200 volts applied to the Ohmjet 
for one to 2 seconds prior to opening of the gas valve, followed by continued power and flow 
for pulse durations of 1.2 seconds. Because of the low gas discharge temperatures attained 
and the limitations on accuracy of thrust measurement, the thrust-time records were essen- 
tially the same for hot flow as for cold flow tests. 

D. Discussion of Ohmjet Testing Under Steady State Operating Conditions - 

In order to obtain data on electrical efficiency and heat transfer coefficients, a series 
of tests were made with the Ohmjet operating at steady state flow conditions. Measurements 
were made of the inlet pressure, P^, the jet discharge total temperature, T^* an d the 
maximum wall temperature, T^^, of the Nichrome heater element. 

The inlet pressure was measured with a Bourdon pressure gage (0 to 100 psig). The 
exit total temperature was measured by impinging the jet into a stagnation cup which was 2 
inches in diameter by 3 inches deep in which a radiation -shielded Chromel-Alumel thermo- 
couple was located. The maximum Nichrome heater temperature for each run was esti- 
mated in the following way. By visual sighting at a slight angle up through the Ohmjet 
nozzle, the heater could be seen and compared in color temperature to that of a 150 watt 
tungsten bulb. The applied voltage to the bulb filament was adjusted to achieve color match 
and the temperature of the filament was determined from a calibration of voltage against 
filament temperature as read by a Pyro -Micro Optical Pyrometer (Model 95). 

The experimental data obtained are shown in Table V-l. Measurements were made 
at applied voltages up to 235 volts with exit CO2 temperatures up to 555° F. and inlet pres- 
sures up to 91 psig. The electrical power input for each point was calculated from the 

measured voltage by using the measured room temperature resistance of 22.55 ohms, and 

o 

assuming a temperature coefficient of resistivity for Nichrome of 0.0004 per C. 
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Figure V-7 shows the experimental data on CO^ discharge total temperature vs. 
electrical power input for a range of inlet pressures from 35 to 105 psia. 

For many of the test conditions, about half the length of the Ohmjet was observed 
to glow cherry red. No thermal insulation was used to reduce radiative or convective heat 
losses. After 13 successful test runs, the Nichrome heater was inadvertently burned out 
at an applied voltage of 180 volts and an inlet pressure of 2 1 psig. The Nichrome failed at 
a point about 1.2 inches upstream from the nozzle. Inspection showed some damage to the 
rear bus -bar and nozzle inlet, probably initiated by fragments or droplets of Nichrome 
which produced an electrical short. 

Since the mass flow rate was not measured directly, the data were reduced using the 
assumptions: (1) that the drop in total pressure through the Ohmjet heat exchanger section 
is negligible and (2) that the mass flow rate corresponds to choked flow (Mach 1) through a 
1/16 inch diameter orifice at a total pressure equal to the measured inlet pressure and a 
total temperature equal to the measured exit total temperature. With these assumptions, 
the mass flow rate (w) for each test point was calculated from the equations: 

w = 0.63 °R 1/2 /sec. for C0 9 (K=1.25) 

A P * 

TH 02 

—5 2 

with P Q2 = P Q1 and A TH = throat area = 2. 13 x 10 ft. 

The power input to the gas was calculated from 

q G =W S <T 02' T 01 ) 

where C = average specific heat of CO„ =0.23 Btu/lb. F. 

P 1 

The electrical efficiency was calculated from 

W (T 02 - V 

= q E 

where q = electrical power input, Btu/sec. 

E 
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The heat transfer coefficient was defined and calculated from 


h = 




a ht aT f 


where the effective heat transfer area is assumed to be one -half the surface area 

of the helical Nichrome heater wire and the film temperature difference (AT_Vis the differ' 
ence between the measured maximum wire temperature (near the nozzle end) and the 
measured gas exit temperature. 

Calculated values of the mass flow rate, electrical efficiency and heat transfer co- 
efficient are tabulated in Table V- 1 . 

The values of electrical efficiency vary from 25 per cent to 55 per cent over the 
range of measurements which does not extend to an inlet pressure of 150 psia and power 
input of 2000 waits for which the unit was originally designed. 

The experimental data show that the heat transfer area should be increased sub- 
stantially in order to obtain the desired gas exit temperature near 1600° F. Based on the 
maximum measured gas discharge temperature of 555° F. in these experiments, the 
specific impulse for expansion to a vacuum with a large area ratio nozzle would be approxi- 
mately 90 seconds. However, it appears to be feasible to attain values of specific impulse 
with OC> 2 of at least 1 10 seconds by proper design of the heat exchanger . 

The heat transfer coefficients were correlated as shown in Figure V-8 by plotting 
against the mass velocity G, which is defined as the mass flow rate divided by the area of 
the pitch circle of the helical Nichrome heater. The data show some scatter and appear to 
vary somewhat more rapidly than the first power of mass velocity within the range of 
measurements. However, the correlation appears quite good considering the accuracy of 
wall temperature measurements (+ 100° F.) and the fact that no corrections have been in- 
troduced for variations in Prandtl number or viscosity of CO~ or for large differences be- 
tween wall and gas temperatures . 
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HEAT TRANSFER COEFFICIENT vs . MASS VELOCITY FOR OHMJET HELICAL HEATER 

FIGURE V-8 
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VI - CONCLUSIONS 


1. For a typical 150 ft. diameter rotating manned space station, it has been found that 
the angular impulse requirements to compensate for gravity gradient torques sub- 
stantially exceed the angular impulse required to precess the spinning station at a 
rate of 360° per year to maintain sun orientation. 

2. The secular gravity gradient torque is the predominant disturbing torque and it is 
periodic, with a period of approximately 54 days, depending on the regression rate 
of the orbit. In a typical worst case, the secular gravity gradient torque, if uncor- 
rected, could misalign a station by 0. 3° per orbit or 4. 8° per day. 

3. The periodic component of gravity gradient torque which averages to zero over each 
orbit, typically produces on angular deviation of the station axis of less than 0. 115° 
which can generally be ignored unless rigid requirements on sun orientation accuracy 
were imposed. 

4 . Corrections for gravity gradient torque can be made by applying a short duration 

(1 to 5 sec . ) thrust impulse either toward the sun or away from the sun at the proper 
point near the station periphery once each station rotation or very N th rotation during 
all or a part of each orbit. The mean angular position of the thrust device during its 
firing, as measured in the plane of rotation^and the net impulse delivered per orbit 
should be controlled and changed gradually over the 54 day cycle of secular gravity 
gradient torque in order to minimize propellant consumption. In this way the correc- 
tive impulse can be applied with the maximum moment arm around the desired axis 
for the corrective torque and the net impulse delivered per orbit can be modulated 
to avoid over -compensation for gravity gradient torque. 

5. Design studies of the Ohm jet electrically-heated pulse rocket system and the Heliojet 
solar -heated pulse rocket system indicate that all the thrust requirements for gravity 
gradient torque compensation, annual precession, and orbit maintenance can be accom- 
plished using as propellant the carbon dioxide generated by the crew if the CO is 
heated to approximately 1600 F . to increase the specific impulse to a value near 120 
seconds . 
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7. Design studies of the Heliojet solar heated pulse rocket system show that two thrust 
units, each having a solar collector diameter of approximately two feet, could accom- 
plish attitude control and orbit maintenance of a 150 ft . diameter space station . The 
system should use a CO recovery and storage system to permit operation of the 
Heliojet every one or two station rotations while it is in sunlight. If the Heliojet were 
designed to operate only for a few minutes per orbit during regeneration of molecular 
sieve beds (which would have to occur in sunlight) the solar collector size required 
would become excessive. 

8. The weight of a complete Heliojet system, including CC> 2 recovery and storage equip- 
ment but excluding the weight of C0 2 consumed, is estimated to be approximately 138 
lbs. as compared to 335 lbs. for an Ohmjet system which uses C0 2 directly from 
molecular sieve beds and 411 lbs. for an Ohmjet system which uses C0 2 compressors 
and storage tanks. Over 60 per cent of an Ohmjet system weight is associated with 

the increased weight of the auxiliarfelectric power system required. For comparison, 
the weight of a chemical thrust system (I = 300 sec.) is estimated to be about 7500 
lbs. for a one year life. 

9. Tests of an Ohmjet thrust unit with pulsed operation at inlet pressures of 125 to 150 
psig have demonstrated thrust levels of 0.31 to 0.375 lbs. with accurate control of 
pulse duration and with repeatability of the impulse delivered per pulse within the 
accuracy of the thrust measuring instrumentation (** 5%) . 

10. Steady state tests of an Ohmjet thrust unit at inlet pressures up to 105 psia and power 
levds up to 1565 watts resulted in gas discharge temperatures up to 555° F . , indicating 
that the heat transfer area would have to be increased to attain an outlet temperature 

of 1600° F. However, steady state electrical efficiencies up to 55 per cent were at- 
tained with no external thermal insulation on the Ohmjet unit. It appears that specific 
impulse values of at least 110 seconds with <X> 2 should be attainable with an improved 
heat exchanger design. 

11. Heat transfer coefficients were measured for the Ohmjet unit which would permit a 
more accurate analysis to optimize the design. 
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VII - RECOMMENDATIONS 


1. In view of the advantages in simplicity and potential reliability of an Ohm jet system 
which uses waste CO directly from a molecular sieve bed, it is recommended that 
Ohmjet thrust units be designed to operate at low inlet pressures near one atmosphere 
and that experiments be conducted to determine the electrical and thermodynamic per- 
formance which is attainable under such conditions. 

2. Additional analytical effort and computer studies should be conducted to determine 
whether there may be optimum orbits (launch date and time of day of launch) which 
would minimize the propellant required for gravity gradient torque compensation on 
space stations for a desired mission duration. 

3. Analytical and design effort should be directed toward specifying design criteria for 

a space station attitude control system which incorporates a gravity gradient computer 
to minimize propellant consumption. 

4. Because of the potential weight saving by using solar power instead of electrical power 
to heat an attitude control propellant, it is recommended that additional design and 
development work be conducted to obtain a more accurate evaluation of the potential 
of a Heliojet attitude control system. Studies indicate that a Heliojet should operate 

at low pressures to simplify or eliminate the problem of CC^ compression and storage. 
In case the compressor can be eliminated, the time for molecular sieve bed regenera- 
tion should be increased so that CO^ can be used directly from sieve beds during a 
large fraction of the time in sunlight, thereby minimizing solar collector size. 
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APPENDIX A 


TRANSIENT HEAT TRANSFER ANALYSIS OF THE HELIOJET 


I. Derivation of Equations - 

The differential equations governing the Heliojet tube wall and gas temperatures 
are derived from appropriate energy balance equations for a differential section. An ap- 
proximate solution is obtained which identifies the significant dimensionless parameters 
and serves as a guide in selecting the preliminary design. 

A. General Assumptions 

1. The tube wall is initially heated externally by solar radiation to a uni- 
form temperature, T , with no gas flow through the tube. 

2 . The duration of gas flow is short so that heat received through the wall 
from additional solar radiation during the gas flow pulse is negligible . 

3. The tube wall is thin so that axial conduction in the wall may be negjdcted 
during gas flow and radial temperature variations across the tube wall 
may be neglected. 

4. The gas flow is essentially incompressible at constant flow rate and 
with a constant heat transfer coefficient. 

5 . The radius of curvature of the tube center line is large compared with 
tube radius so the flow may be analyzed on the basis of one dimensional 
flow in a straight tube . 

6. The flow is quasi-static; hence the gas temperature variation with dis- 
tance along the tube may be found as a steady-state distribution correp 
sponding to the instantaneous wall temperature . 


B. Nomenclature 

2 

A = tube cross-sectional area 7D 

4 

C C = specific heats for gas and wall 
D = tube diameter 
G = gas mass velocity (=pV) 
h * convection heat transfer coefficient 


BTU/lb° F. 
ft 

lb/ft 2 sec 
BTU/ft 2 sec° F 
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K i =4 (sr^)(i5) =4S, (i) 

' Pg 


ft 


-1 


= K 1 L 


K„ = 


6p C 
K w pw 


= St 


fe) te) 

pw 


K, = 


-K ' 
1 - e 2 


K r 


L = tube length 

_ h 

St = Stanton Number = — — ^r 

Pg 

t = time 

T t , T , = gas temperature at any point, at exit, 
g’ ge' gi 

at inlet 

T' = gas temperature averaged over tube length 

t-v X t = wall temperature at any point, averaged 
w’ wa, wo 

over tube length, initial uniform tem- 


perature 

V = gas velocity 
w = weight flow rate of gas 
x = distance along tube from inlet 
6 = tube wall thickness 


p > P 
M g w 


= weight density of gas, tube wall material 


sec 


-1 


sec 


ft 


-1 


sec 


F 

ft/sec 

lb/sec 

ft 

ft 

lb/ft 3 


C. Heat Balance Equations 

In view of assumptions 2, 3, and 5 the heat exchanger tube is treated as 
an externally insulated straight tube with wall temperature as a function of x and t as 
shown in Figure A.-1. 
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tube geometry for heat transfer analysis 

FIGURE A- 1 

A unit mass of gas moving a distance dx in unit time acquires from the wall thermal 
energy equal to C dT ; hence for w = p VA = GA the heat balance equation is: 

G-j- D 2 C dT = hnDdx (T - T ) 

4 pg g w g 

HT _ 4 hdx _ T v (1) 

dT g C DG ( w g 
6 Pg 


Dividing by d t and using: 


dT *T 

S' r C. 


aT *t 


+ w -7T = 


we get: 1 ^ T g 

V *t + 


— — — (T - T J 
C DG v w g 
Pg 


Consider next the heat loss from the wall which equals the rate of decrease of 

stored energy in the wall, i.e.: -v-p 


^T 

p C rrD dx = -hnDdx (T - T ) 

/ W pW Q t w * 


W g' 


(T - T 

W £ 
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ditions: 


Equations (2) and (3) must be solved simultaneously subject to the boundary con 


at t - o: T (x, o) = T = const, 
w wo 

at t = o: T (o,t) = T = const. 

S g 1 


In view of the quasi -static assumption, 6: 




-L Zjs. « 

V Jt 


£T 


JS. 


(4) 


and the two equations to be solved may be written as: 

* T g 

sr - K i y - y 


a t 


= ~ K s (t w - y 


(5) 


( 6 ) 


The Stanton number, which is defined to befh/C d is assumed to be given by 
the empirical correlation: Pg 


St = 


C G 

Pg 


= 0.023 Re' 0 * 2 Pr ' 0 * 6 


where Re = Reynolds number 
Pr = Prandtl number 

This equation is used to compute the heat transfer coefficient, h, for flow of CO inside 

tubes. The effect of the film temperature difference on the heat transfer coefficient has 
not been included for this first order analysis. 

II- Solution of Equations - 

An upper limit for exit temperature, T ge at t * o may be obtained by using the 
initial wall temperature T wq in place of T w in (5). The upper limit for T may then be 
found by solving (6) with T g = Tge (0) and using this value of T w « in (5). *ki result of 


such an approach is: 


AT 


T (t) 


max 


- gl 


ge max 


T -T , 
wo gl 


= !— e "*2 


[l+(l- 


e K 2) (1 -e~ K 3 t 


)] (7) 


where: K 2 = K^L 


77 



This solution *ives a reasonable value of T (0), assuming that the starting transient is 
short and wall temperature does not drop appreciably during this transient. However, 


equation (7) yields too high a gas exit temperature at t = cn. 

A probable lower limit may be obtained by solving (6) with 
temperature instead of the exit temperature as above. For this case 


T (t) 

ge 


AT . (t) 
ge min 


mm 


gl 


= (1 — e ~^2) 


e-V 


wo 


gl 


= T , , the inlet 
gl 

the result is: 




Tliis solution starts (t = o) at the same point as the maximum limit case, but falls off 
too rapidly to the proper final value of T gg (oo) = T^. 

A much better estimate of T .(t) lying between the two limits above may be made 
by average values for the temperature in solving (5) and (6). This solution is worked 

out in detail below. L 

In equation (5) use = — J" T^dx, i.e. the heat flux is based upon an 

average wall temperature instead of the extreme case (T w = T wq ) as in the previous 
solutions . 


^ T 

. = K, (T -T ) 

Jx 1 wa g 

Equation (10) has the solution: 


L T w*>- T g J e ” lc ‘ x 


where: at x = o: T (0, t) = T . * inlet temperature 

g g 1 

at x = oo :T (oo, t) = T = average wall temperature, 
g w® 

Equation (6^ if averaged over L is: 

* T 

9 wo 

■ ■■■- — m 

From (11) we have: 


K 0 (T - T ) 
3 wa ga 


( 10 ) 


( 11 ) 


( 12 ) 


T (t) = — 
ga' L 


y \ <*. t) dx = Twa « - i <1 - .““l) fT wa (t) 

o © 2 




(13) 
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Substituting into (12) we get: 


(-■*) v-S (-■*■) 


with solution: 


T (t) = T , + (t -T .) e 
wa gl \ wo gl/ 


-K 4 * 


g 1 ' 


K 3 /, 

where: K 4 =— (l-e ) 


at t = o : T (0 ) = T = initial wall temperature 

wa wo 

at t = oo : T (oo) = T , = inlet gas temperature 
wa gl 


We now substitute (15) into (11) to get: 

T (x, t) = T . + - (t -T . ) ( l-e 1 ) 
e ffl v. wo irl • ' ' 


gi 


-V 


e 


(14) 


(15) 

(16) 


(17) 


— K,x 

m Ui- 

g ' ' gl ^ wo gl 

which checks all the initial and final conditions. The exit temperature is accordingly: 



T (t) 
ff : 

T 

wo 




(18) 


Repeated iterations do not change (17); hence, the answer is the best possible based on 
the assumption that the wall is isothermal at each instant of time. Equation (18) is 
plotted in Figure A -2. The average wall temperature for equation (15) is shown in Figure 
A-3. A comparison of all solutions for the special case =3.0 is given in Figure A-4. 

Exact analytical solutions for equations (5) and (6) which have been found by the 
use of Laplace transforms are too unwieldy for easy use. 
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